Episodic synchronization in dynamically driven neurons 
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We examine the response of type II excitable neurons to trains of synaptic pulses, as a function 
of the pulse frequency and amplitude. We show that the resonant behavior characteristic of type II 
excitability, already described for harmonic inputs, is also present for pulsed inputs. With this in 
mind, we study the response of neurons to pulsed input trains whose frequency varies continuously in 
time, and observe that the receiving neuron synchronizes episodically to the input pulses, whenever 
the pulse frequency lies within the neuron's locking range. We propose this behavior as a mech- 
anism of rate-code detection in neuronal populations. The results are obtained both in numerical 
simulations of the Morris-Lecar model and in an electronic implementation of the FitzHugh-Nagumo 
system, evidencing the robustness of the phenomenon. 
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I. INTRODUCTION 

Neurons exhibit all-or-none responses to external input 
signals. The main function of this thresholding behavior 
is to process information in a way that is efficient and 
robust to noise p). Input signals received by most non- 
sensory neurons take the form of pulse trains, coming 
from the spiking activity of neighboring neurons. There- 
fore, in order to understand the mechanisms of informa- 
tion processing in neural systems, it is very important 
to characterize in detail the response of neurons to pulse 
trains. Furthermore, realistic pulse trains are intrinsi- 
cally dynamical, with an instantaneous firing frequency 
that varies continuously in time. It is therefore neces- 
sary to assess the influence of this non-stationarity in the 
neuronal response. This paper addresses these questions. 

Most studies of driven neurons have been restricted so 
far to harmonic driving signals 0, IS E3 IE @ • Many of 
these works have shown that for certain types of neurons, 
i.e. those exhibiting what is called type II excitability, a 
resonant behavior arises with respect to the external driv- 
ing frequency. Excitability in those neurons is usually as- 
sociated with an Andronov-Hopf bifurcation, which leads 
to the existence of subthreshold oscillations in the ex- 
citable regime. When the frequency of these oscillations 
equals that of the harmonic driving, a resonance arises. 

It is to be expected that a similar resonant behavior ex- 
ists for pulsed inputs. In that case, the same pulse train 
impinging on two different neurons could elicit a response 
on only one of them, i.e. on the one that is tuned to res- 
onate with the incoming pulse frequency. This behavior 
has indeed been observed experimentally in a rat's neo- 



cortical pyramidal neuron that innervates another pyra- 
midal neuron and an interneuron; a bursting input from 
the innervating neuron produced an action potential in 
the interneuron but not in the second pyramidal neuron 
• This behavior was interpreted in terms of a differen- 
tial frequency-dependent facilitation and depression, re- 
spectively, and as such was studied by Izhikevich and 
co-workers Here we propose a simpler mechanism 

for this phenomenon, relying only on the resonant be- 
havior of the processing neuron. This mechanism could 
provide a means for distinguishing between firing rates, 
whose controlled variation lies at the heart of the rate 
coding approach to information processing by neurons. 

Our results show that type II excitable neurons ex- 
hibit a resonant response with respect to the frequency 
of input pulse trains. This behavior leads to episodic 
synchronization between the neuron's output and an in- 
put with dynamically varying firing rate. Episodic syn- 
chronization has previously been reported in coupled 
lasers with intrinsic dynamics |10|. Here we extend that 
property to externally driven excitable systems. Two 
types of systems have been investigated: a Morris-Lecar 
model (Sec. [HJ and an electronic implementation of the 
FitzHugh-Nagumo model (Sec. IIIII) . 



II. MORRIS-LECAR MODEL 

A. Model description 

We consider neurons whose dynamical behavior is de- 
scribed by the Morris-Lecar model |Tl|. 
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where V and W represent the membrane potential and 
the fraction of open potasium channels, respectively. C m 
is the membrane capacitance per unit area and <j) is the 
decay rate of W . The neuron is affected by several cur- 
rents, including an external current / a pp, a synaptic cur- 
rent I S ym and an ionic current given by 

Iion = g C aM 00 (V)(V-V$ a ) + 

gKW(V-V°)+g L (V-V°). (3) 

In this expression, g a (a — Ca, K, L) are the conduc- 
tances and V a ° the resting potentials of the calcium, 
potassium and leaking channels, respectively. We define 
the following functions of the membrane potential: 
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where Vmi, Vmi, Vwi and Vw2 are constants to be spec- 
ified later. The last term in Eq. is a white Gaussian 
noise term of zero mean and amplitude D. 

In the absence of noise, an isolated Morris-Lecar neu- 
ron shows a bifurcation to a limit cycle for increasing 
applied current I app Depending on the parameters, 
this bifurcation can be of the saddle-node or the subcrit- 
ical Hopf types, corresponding to either type I or type 
II excitability, respectively. The specific values of the 
parameters used are shown in table [I] 12]. For these 
parameters, the threshold values of the applied current 
under constant stimulation are 39.7 mA for type I and 
46.8 mA for type II. 

In this paper we analyze the behavior of a neuron 
driven by a synaptic current. To that end, we use the 
simplified model of chemical synapse proposed in Ref. 
|13| , according to which the synaptic current is given by 



g syn r(t)(V - E„), 



(7) 



where g syn is the conductance of the synaptic channel, 
r(t) represents the fraction of bound receptors, and E s is 
a parameter whose value determines the type of synapse: 
if E s is larger than the rest potential the synapse is exci- 
tatory, if smaller it is inhibitory; here we consider an ex- 
citatory synapse with E s = mV. The fraction of bound 
receptors, r(t), follows the equation 



^ = a[T](l-r)-(3r, 



(8) 



where [T] = T max 8(T + r syn - t)6(t - T ) is the con- 
centration of neurotransmitter released into the synaptic 
cleft by the presynaptic neuron, whose dynamics is also 
given by Eqs. CQ)-I0 with no synaptic input, a and /3 
are rise and decay time constants, respectively, and Tq 
is the time at which the presynaptic neuron fires, which 



happens whenever the presynaptic membrane potential 
exceeds a predetermined threshold value, in our case cho- 
sen to be 10 mV. This thresholding mechanism lies at 
the origin of the nonlinear character of the synaptic cou- 
pling. The time during which the synaptic connection is 
active is roughly given by T SV n . The values of the cou- 
pling parameters that we use [13j are specified in Table 
IB The equations were integrated using the Heun method 
(l4| . which is a second order Runge-Kutta algorithm for 
stochastic equations. 



Parameter 


Morris-Lecar Til (TI) 


c m 


5/iF/cm 2 




8/iS/cm 2 


9l 


2 /iS/cm 2 


9C a 


4.0^S/cm 2 


V K 


-80 mV 


Vl 


-60 mV 


V Ca 


120 mV 


Vmi 


-1.2 mV 


Vm2 


18 mV 


Vwi 


2mV (12 mV) 


V\V2 


17.4 mV 




1/15S- 1 


Parameter 


Synapse 




a 


2.0ms~ 1 mM" 1 




1.0ms" 1 


Tmax 


1.0 mM 


ffsyn 


(specified in each case) 


T syll 


1.5 ms 


E s 


OmV 



TABLE I: Parameters values of the Morris- 
Lecar and synapse models used in this work. 



Response diagram of a periodically driven 
driven neuron 



First we analyze how a Morris-Lecar neuron responds 
to periodic inputs of varying frequencies. Specifically, we 
ask how large the signal needs to be in order to elicit 
spikes in the receiving neuron. It is also important to 
characterize the frequency of spiking in terms of fre- 
quency of the input. As mentioned in the introduction, 
this question has already been addressed, in the case of 
harmonic inputs, for different neuronal models, including 
the Morris-Lecar model |lf| ■ We will now compare these 
results with those obtained for a pulsed input. For the 
Morris-Lecar model, one can expect a completely differ- 
ent behavior between the type I and the type II cases, 
given that the bifurcation to a limit cycle is a saddle- 
node bifurcation in the former case and a Hopf bifurca- 
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tion (with the well known eigenfrequency associated to 
the spiral fixed point) in the latter. 

We first consider an isolated neuron without synaptic 
inputs (J S yn — 0), but subject to a harmonic modulation 
of the applied current 7 app with the form, 



7 + Acos(27r/ in t). 



(9) 



In order to quantify the response of the neuron to this 
harmonic input, we plot in Fig.^(in color scale) the ratio 
between the output and the input frequencies, / ou t /fin, 
as a function of the amplitude A and frequency f m of the 
applied current ©. The figure compares the resulting 
response diagrams of type I and type II neurons. To 
obtain these plots, A is varied for fixed /, while using 
as initial condition for a given A the final state of the 
previous A value. In the upper panels A increases, thus 
showing the stability of the rest state, while in the lower 
panels A decreases, this indicating the stability of the 
limit cycle. The figure shows that type II neurons have a 
region of bistability, where the fixed point and the limit 
cycle coexist. In contrast, for type I neurons the two 
plots are basically the same, indicating and absence of 
bistability. 




FIG. 1: Response diagram of Morris Lecar neurons for a 
harmonic input: /out //in is plotted (color-coded) as a function 
of the amplitude and frequency of the applied current. Left 
plots: type I neuron (with 7o = 39 mA), right plots: type 
II neuron (with 7o = 46 mA). The response is measured for 
increasing A in the upper plots, and for decreasing A in the 
lower plots. 

There is another qualitative difference between types 
I and II that can be observed in Fig. ^ In the type 
I neuron, the critical modulation amplitude for spiking 
increases monotonically with the frequency of the stimu- 
lus. On the other hand, in the type II neuron the critical 
amplitude exhibits a minimum for a given nonzero fre- 
quency, in our case around 20 Hz. This behavior can be 
understood as resulting from the subthreshold damped 
oscillations characteristic of type II excitability . 



We now characterize the response of a neuron to an in- 
put train of periodic synaptic pulses of varying frequen- 
cies and amplitudes. To that end, we drive the neuron 
with a synaptic current with the form 



I syn = Ar{t,f ijl )(V(t)-E s ). 



(10) 



In this model, the frequency is given by the dynamics of 
r(t, /in), described in Eq. (JSJ, assuming that the presy- 
naptic firings occur periodically with frequency /; n . In 
this way, we can quantify the response of the neuron in 
terms of the efficiency in responding to a periodic synap- 
tic input with a given frequency and amplitude, as shown 
in the harmonic case. Figure [21 shows the corresponding 
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FIG. 2: Response diagrams of Morris Lecar neurons to a 
periodic synaptic (pulsed) input: /out //in is plotted (in color 
scale) as a function of the amplitude and the frequency of 
the synaptic current. Left plots: type I neuron (with 7 apP = 
39 mA), right plots: type II neuron (with 7 app = 46 mA). The 
response is measured for increasing A in the upper plots, and 
for decreasing A in the lower plots. 

response diagrams, i.e. /out //in as a function of A and 
/in, for both excitability types and for increasing (top) 
and decreasing (bottom) A. The behavior shows features 
common to the harmonic case, such as the existence of 
the same resonant frequency in type II for both kinds 
of inputs. But there are also very interesting differences 
between them, specially in the high and zero frequency 
limits. 

The main difference between the harmonic and pulsed 
input cases is the approach to the DC threshold current 
(39.7 mA for type I and 46.8 mA for type II). While in 
the former case this happens for frequencies approaching 
zero (type I) or resonance (type II), for pulsed inputs it 
happens for high frequencies, i.e. when the signal period 
is of the order of the pulse width. This is the reason for 
the appearance of a spiking region at high frequencies 
for pulsed inputs, which is absent in the harmonic case 
|l6j| . Also, in the low- frequency limit one can observe, for 
pulsed inputs, a constant value of the critical amplitude. 
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This is related with the fact that, when the input pe- 
riod is high enough with respect to the pulse width, the 
system response is essentially independent of the period. 



C. The dynamical case: Variation of the input 
frequency 

In the previous section, we have characterized the be- 
havior of a neuron subject to a pulsed synaptic current 
of fixed frequency. Our results show that type II neurons 
exhibit a resonant behavior, defined by the existence of 
an optimal frequency for which the critical amplitude for 
spiking is minimal. The question now is, what happens if 
the frequency of the input train varies dynamically, which 
is a more realistic situation for a non-sensory neuron. 

To answer this question, we made simulations with two 
Morris-Lecar neurons coupled unidircctionally through a 
chemical synapse. The input neuron operates in the limit 
cycle regime and is considered to be type I, so that we 
can control its spiking frequency by varying its applied 
current / app 0. This neuron is synaptically coupled to 
a type II neuron operating in an excitable regime, with 
a coupling strength <7 syn such that the receiving neuron 
only fires in a given range of frequencies (i.e. the cou- 
pling is such that the amplitude of the input pulses lies 
below the critical amplitude at zero frequency but above 
its minimum at resonance; this corresponds e.g. to a hor- 
izontal line at around 0.4 mV in the right plots of Fig.EJ. 

Figurc|3shows what happens when the firing frequency 
of the input neuron first increases and then decreases in 
the range 13-28 Hz. The plot compares the instantaneous 
firing frequency of both neurons, relating them with the 
boundaries of the locking range of the second neuron, in- 
dicated by horizontal dashed lines; the 1:1 locking region 
is specifically shown. It can be seen that as the input fre- 
quency increases (first half of the plot) , the receiving neu- 
ron starts spiking with approximately 1:1 frequency ratio 
when the input frequency falls within the corresponding 
locking range, the ratio decreasing when the input fre- 
quency exceeds ^20 Hz. Spiking persists while the input 
frequency remains in the wider (not 1:1) locking range, 
and is maintained even for a while after the input finally 
exits the locking region. A similar behavior is observed 
for decreasing frequencies, but the "inertia" observed at 
the exit of the locking region is larger than for increasing 
frequencies. The time series of the receiving neuron is 
shown in the lower plot; the episodes of synchronization 
with the input signal are clearly observed. 

In order to understand the dynamic driving effects re- 
ported above, and particularly the locking asymmetry 
observed between an increase and a decrease in the in- 
put firing frequency, we now study the neuron response to 
a controlled variation of the frequency for different vari- 
ation rates. To that end, we consider a synaptic input 
whose frequency is uniformly changing from one cycle to 
the next at a rate /' = A// At, and measure the response 
of the neuron in terms of its instantaneous output fre- 
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FIG. 3: Upper plot: instantaneous frequency versus time of 
the input neuron (black stars) and of the receiving neuron 
(red triangles). The three horizontal dashed lines indicate 
the boundaries of the locking region for the selected synap- 
tic strength (g Byn = 0.43 nS); the 1:1 region is specifically 
shown, and labeled at the right. Lower plot: time series of 
the receiving neuron. 

quency / ou t- Figure 01 shows the response diagram for a 
fixed synaptic strength and different frequency variation 
rates, for both increasing and decreasing input frequen- 
cies /i n . The figure shows that the slower the variation 




FIG. 4: Response of a type II neuron to a controlled variation 
of the input frequency, for two different rates of this variation, 
/' = Af /At, for g syn = 0.38 nS, compared with the adiabatic 
response (which is a horizontal cut of the response diagram 
at A = 0.38 nS). 

rate, the closer the response is to an adiabatic passage, 
as expected. Additionally, the results indicate that the 
persistence of the output neuron in the firing state (even 
when the input signal has left the locking region) is much 
larger when the frequency decreases than when it in- 



creases. This is consistent with the asymmetric response 
exhibited in Fig. [3J and can be expected to arise from 
the asymmetric shape of the response function / ut//in ; 
which is equal to 1 for small frequencies and moderately 
smaller than 1 for large frequencies. Evidently the neu- 
ron prefers to respond in a 1:1 regime, which produces a 
larger persistence for decreasing frequencies. 

To further quantify the approach to the adiabatic re- 
sponse in terms of the rate of change in the input fre- 
quency, we can define a distance T> to this adiabatic re- 
sponse as the absolute value of the difference between 
the area of the neuron's response diagram at a given rate 
/', as plotted in Fig. and the area of the adiabatic 
response: 



V = 



fn 



fin 
/out 



out / adiab 



dfu 



(11) 



This measure is plotted in the inset of Fig. 0] as a func- 
tion of the rate of change in the input frequency. The 
plot shows that the distance increases as the frequency 
changes more rapidly, as expected. 

The dynamical response described in the previous 
paragraphs leads to episodic synchronization when the 
input pulse train exhibits a varying firing rate. This sit- 
uation is shown in Fig. in which an input pulse train 
whose firing rate takes the form, by way of example, of an 
Ornstein-Uhlenbeck noise with amplitude A ou = 0.1 mA 
and correlation time r ou = 1 s in the ~ 20 — 36 Hz range. 
The response of the second neuron for g syn = 0.38 nS is 
displayed in the bottom plot, and exhibits clear episodes 
of synchronization with the input signal, whenever the 
firing rate of the later falls (approximately) within the 
locking range of the neuron for the coupling strength cho- 
sen (represented by horizontal dashed lines in the figure) . 
In that way the receiving neuron acts as a bandpass filter 
for input pulse trains. 



III. FITZHUGH-NAGUMO CIRCUIT 
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FIG. 5: Upper panel: Instantaneous frequency of the synaptic 
input (black stars) and of the output neuron (red triangles) 
as a function of time. The horizontal dashed lines delimit the 
region of locking according to the adiabatic calculations of 
Fig-IUfor the synaptic strength used (<? S yn = 0.38 nS). Lower 
panel: time series of the input signal. 
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FIG. 6: Response diagram of an electronic FHN neuron for 
a periodic input of pulses with 10 ms width. The limit of the 
region of 1:1 resonance is marked with a dashed line. 



In order to show that the behavior reported in the pre- 
vious Section is generic and robust, we have reproduced 
the results with an electronic neuron, specifically with 
an electronic implementation of the FitzHugh-Nagumo 
(FHN) model |l7|. The circuit has been previously de- 
scribed in |18| , where synchronization between two FHN 
neurons was studied. A detailed description of the cir- 
cuit can be found at 

Hi- 

In our particular setup, a FHN 
neuron is excited by a pulsed input of variable frequency. 

Following the procedure of Sec. [DJ we first determine 
the response of the electronic neuron to a train of periodic 
pulsed inputs of fixed frequency. The pulses have the 
form of square pulses of 10 ms width. Figure shows the 
corresponding response diagram, obtained by increasing 
the amplitude of the input pulses until the neuron starts 
firing. Similar results (not shown here) are obtained with 
pulses of different width. At first glance, we can observe a 



resonance minimum around /;„ = 15 Hz, which confirms 
that the FHN neuron is of type 11. Two local minima are 
also observed around /;„ = 7.5 Hz and /j n = 31 Hz. It 
is worth noting that despite the spike threshold is low, 
moderately large values of the input voltage are required 
to induce spiking at the input frequency (see region 1:1 
in Fig. El). 

We have thus a type II electronic neuron that exhibits 
a resonance at a frequency close to 15 Hz. Following 
again the approach of Sec.[H] we now subject the circuit 
to pulse trains with time-varying frequency. Specifically, 
the pulse frequency is made to depend linearly with time 
(with both positive and negative slope). Similar results 
(not shown here) are obtained with sinusoidal variations. 
The neuron response to this dynamical input is shown 
in Figs. 0and|Hl where two different input voltages, cor- 
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responding to the values denoted as A and B in Fig. 
have been applied. 
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FIG. 7: Response of a FHN electronic neuron to a pulsed 
input of variable frequency. The input voltage corresponds 
to the value marked by A in Fig. |S| Plot (a) shows the in- 
stantaneous frequency of the input pulses (solid line) and of 
the neuron's output (stars). The shaded region corresponds 
to the frequency ranges for which locking should occur. Plot 
(b) displays the time evolution of the membrane voltage of 
the electronic neuron, which corresponds with the voltage U2 
at condenser CI, in the circuit given in |Tfll |. 




time (s) 

FIG. 8: Response of a FHN electronic neuron to a pulsed 
input of variable frequency. The input voltage corresponds to 
the value marked by B in Fig. |S| Figure layout is as described 
in the caption of Fig. Q 

In the case of Fig. [7| the input signal scans the region 
marked with A in Fig. The upper plot shows the in- 
stantaneous frequency of the input train (solid line) and 
of the FHN neuron (stars). Figure|7fb) plots the neuron's 
output. The results show that the neuron pulses when 
the input frequency lies within the resonance regions 
given in Fig. EI and highlighted in gray in Fig.0[a). This 



behavior is an agreement with the observations made in 
the Morris-Lecar model. The fact that no inertia effects 
are seen when the input frequency sweeps past the reso- 
nance region is due to the frequency variation rate being 
very slow with respect to the characteristic time scales of 
the system (adiabatic limit). 

Figure [8] shows the system's behavior for a different 
value of the input voltage (marked as B in Fig. |SJ, for 
which the input frequency encounters three resonance re- 
gions as it varies. Accordingly, the electronic neuron fires 
whenever the input frequency lies inside any these re- 
gions, exhibiting clear episodes of synchronization. In 
other words, the neuron acts as a band-pass filter, with a 
frequency range that depends on the input voltage level 
according to its response diagram. 

IV. DISCUSSION 

Neurons are information-processing devices. The na- 
ture of coding in neuronal systems is still an open ques- 
tion. One of the most favored views in the field is that 
of rate coding, where the intensity of a signal is encoded 
in the firing rate. Neuronal systems must therefore be 
able to distinguish between firing rates. We have pro- 
posed a way to accomplish that, through the resonant 
behavior exhibited by type II neurons. A population of 
neurons with different tuning characteristics, and there- 
fore distinct locking ranges, should be able to distinguish 
between different incoming pulse frequencies by activat- 
ing selectively different subpopulations that respond se- 
lectively to different frequencies. 

We have systematically analyzed the response of type 
I and II neurons to pulsed driving, compared it with the 
standard case of sinusoidal driving, and observed the res- 
onant behavior of type II neurons. This phenomenology 
leads to episodic synchronization between the input and 
the output of the neuron, when the input consists of a 
train of pulses with dynamically varying frequency. The 
phenomenon has been reported both in numerical simu- 
lations of the Morris-Lecar model, and in an experimen- 
tal implementation of the FitzHugh-Nagumo circuit. We 
expect this type of behavior to underlie information pro- 
cessing in rate coding and decoding neuronal populations 
in more complex brain networks. 
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